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ABSTRACT 


To add genetic information to the conservation efforts on grey goral (Naemor- 
hedus caudatus) in Korea, we investigated the pattern of mitochondrial cytochrome 
b gene sequence (606 bp) of six specimens from two localities in Korea. The 
corresponding sequences of N. caudatus in China obtained from GenBank were 
also used. The nucleotide Tamura-Nei distances between each of four haplotypes of 
N. caudatus in Korea and the haplotype of N. caudatus in China varied from 
0.0650 to 0.0803: N. caudatus revealed high level of sequence diversity in 
Bovidae. In N. caudatus in Korea, the distances among three haplotypes at Yanggu 
were 0.0151 to 0.0185, and it suggests that the genetic diversity of Yanggu 
population was decreased in low level. Moreover, the distances between each of 
three haplotypes at Yanggu and one haplotype at Samcheok were 0.0343 to 
0.0489. It indicates that habitat isolation caused the continuous increase of genetic 
distance with geographic distance in N. caudatus, and various conservation plans 
for mitigating the loss of genetic diversity in Korea have to be in immediate action. 
To clarify the taxonomic status of N. caudatus, the sequence (276 bp) of N. goral 
available from GenBank were also utilized, and N. goral was not distinct from N. 
caudatus. It suggests that they may be conspecific, but further analyses with 


additional specimens of two species are necessary. 
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INTRODUCTION 


Genetic studies in endangered species have become widespread in the past decade, and 
conservation applications have transformed population genetics from academic discipline to the 
forefront of global or local management decisions (Hedrick, 2001). Furthermore, mitochondrial 
DNA (mtDNA) is a highly sensitive genetic marker suitable for studies of closely related taxa or 
populations of a variety of species because of its fast rate of evolution and characteristic maternal 
inheritance (Wilson et al., 1985). Mitochondrial cytochrome b gene analysis has been used to infer 
molecular phylogeny in Bovidae (Chikuni et al., 1995; Groves and Shields, 1996; Birungi and 
Arctander, 2001), and mtDNA has been widely used in genetic studies of population differentiation 
(Avise, 1986). Mitochondrial DNA control region sequences from six Kenyan Grant’s gazelle 
(Gazelle granti) populations were highly divergent (about 14%) among populations (Arctander et 
al., 1996). 

Grey goral (Naemorhedus caudatus) distributes in Soviet Far East, eastern China, eastern 
Burma, and western Thailand (Grubb, 1993). In Korea, this species has been protected as one of 
natural monuments, and the Ministry of Environment classifies it as endangered. However, the 
taxonomy of genus Naemorhedus is still in confusion. Nowak (1999) recognized three species of 
gorals (N. baileyi, red goral from Tibet to Burma; N. caudatus, grey goral from southeastern 
Siberia to Burma and Thailand, and N. goral, Himalayan goral from northern Pakistan to Assam) 
and three species of serows (Capricornis sumatraensis, Sumatran serow; C. swinhoei, Formosan 
serow; C. crispus, Japanese serow). These six species were classified into a single genus 
Naemorhedus (Grubb, 1993), but cytochrome b sequence analyses supported the classification of 
two separate genus (Groves and Shields, 1996). Moreover, N. caudatus and N. baileyi were 
subspecies of N. goral (Ellerman and Morrison-Scott, 1951). 

To add genetic information to the conservation efforts of grey goral (N. caudatus) in Korea and 
to clarify the taxonomic status of this species, we investigated the pattern of mitochondrial 


cytochrome b gene sequence of six specimens from two localities in Korea. 


MATERIALS AND METHODS 


Six specimens of grey gorals (N. caudatus) were collected from two localities (2, Samcheok; 4, 
Yanggu) in Korea, and muscles and scat were kept in a deep-freezer or 9596 ethyl alcohol. Total 
cellular DNA was extracted from muscles and scat (Hillis et al., 1996). The cytochrome b gene was 
PCR-amplified using two sets of primers [L14724, H15149; L15513, H15915] designed by Irwin 
et al. (1991). PCR thermal cycle profiles were as follows: 94°C for 5 min, 94°C for 1 min, 57°C for 
1 min, 72°C for 1 min (32 cycles), 72°C for 5 min. For sequencing, the purified PCR products were 
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analyzed with an automated DNA sequencer (model 377, Perkin Elmer) at Intron Co. (Seoul, 
Korea). Two parts of sequences (276 bp and 330 bp) in the mitochondrial cytochrome b gene of N. 
caudatus were aligned by eye and no ambiguity in the alignment was observed. The corresponding 
two parts (sequence nos. 156-431 and 811-1140) of sequences of N. caudatus in China obtained 
from GenBank (accession number U17861) were used. In N. goral, however, only the first part of 
sequences available from GenBank (accession number D32196) was utilized. The two parts of 
sequence of Saiga tatarica from GenBank (accession number AF064487) were used as an 
outgroup. Tamura-Nei distances (Tamura and Nei, 1993) were calculated and phylogenetic trees 
were constructed by neighbor-joining and maximum parsimony methods with 500 bootstrap 
replications by the program MEGA (version 1.02). 


RESULTS 


The haplotypes I to IV (276 plus 330 bp) of mtDNA cytochrome b gene in N. caudatus in Korea 
are shown in Fig. 1: haplotype V (276 plus 330 bp) of N. caudatus in China and haplotype VI 
(276 bp) of N. goral from GenBank are also shown. The number of specimens and locality of four 
haplotypes in N. caudatus from Korea are as follows: I, 2 from Samcheok; II, 1 from Yanggu; III, 
2 from Yanggu; IV, 1 from Yanggu. 

Nucleotide Tamura-Nei distances among five haplotypes of mtDNA cytochrome b gene of N. 
caudatus in Korea and China are given in Table 1. In N. caudatus in Korea, the nucleotide 
Tamura-Nei distances among three haplotypes (II to IV) at Yanggu were 0.0151 to 0.0185, and 
the distances between each of three haplotypes (II to IV) at Yanggu and one haplotype (I) at 
Samcheok were 0.0343 to 0.0489. Moreover, the distances between each of four haplotypes (I to 
IV) of N. caudatus in Korea and the haplotype (V) of N. caudatus in China varied from 0.0650 to 
0.0803. 

Phylogenetic trees constructed with five haplotypes of mtDNA cytochrome b gene (276 plus 330 
bp) of N. caudatus in Korea and China are shown in Fig. 2 (A, neighbor-joining tree; B, maximum 
parsimony tree). The sequence (S) of Saiga tatarica from GenBank was used as an outgroup. Three 
haplotypes (Il to IV) at Yanggu differed from haplotype I at Samcheok, and four haplotypes (I to IV) 
in Korea was different from haplotype V in China. In N. caudatus, genetic diversity appeared to be 
decreased at Yanggu population, and genetic distances were increased with geographic distance. 

Moreover, the trees with six haplotypes of mtDNA cytochrome b gene (276 bp) of two species in 
Naemorhedus are shown in Fig. 3. Haplotype VI of N. goral formed a subgroup with two 
haplotypes of N. caudatus (I, Samcheok, Korea, and V, China). In cytochrome b gene, N. goral 
was not distinct from N. caudatus. 


DISCUSSION 


In Bovidae, mtDNA cytochrome b gene exhibited considerable infraspecific variation, i.e., 12- 
1496 in Grant's gazelle, Gazella granti (Arctander et al., 1995) and 0.5-1.196 in lechwe, Kobus 
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I AATACACTAT ACATCTGACA CAACAACAGC ix ed GTAACACACA TCTGCCGAGA 215 
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I AAACAACGAA GCATGATATT CCGACCAATC AGCCAATGCT TATTCTGAAC TCTAGTAGCA 990 


I GGACAACTGG CATCCATCAT ATATTTCTTC ATCATCCTGG TACTAATACC AGTAGCTAGC 1110 


Fig. 1. Nucleotide sequence alignments (276 plus 330 bp) of the mitochondrial cytochrome b gene in two 
species of Naemorhedus. In N. caudatus, five haplotypes of all 606 bp were shown (haplotypes I to IV, Korea, 
this study; haplotype V, China, GenBank). The number of specimens and locality of four haplotypes from 
Korea are as follows: I, 2 from Samcheok; Il, 1 from Yanggu; M, 2 from Yanggu; IV, 1 from Yanggu. In N. 
goral, haplotype VI of 276 bp available from GenBank was shown. Arabian numerals indicate the sequence 


numbers of cytochrome b gene of N. caudatus in GenBank. 
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Table 1. Tamura-Nei distances among five cytochrome b haplotypes in Naemorhedus caudatus in Korea 


and China. Haplotypes are labeled as in Fig. 1. 


Types I II Il IV 
II 0.0343 
Il 0.0489 0.0185 
IV 0.0486 0.0168 0.0151 
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Fig. 2. Phylogenetic trees with five haplotypes (606 bp) of Naemorhedus caudatus from Korea and China. 
Trees were constructed by neighbor-joining (A) and maximum parsimony (B) with 500 bootstrapped 
replications. Tamura-Nei distances presented in Table 1 were used. Haplotypes are labeled as in Fig. 1, and 


the sequence (S) of Saiga tatarica was used as an outgroup. 


leche, and 0.696 in waterbuck, K. ellipsiprymnus (Birungi and Arctander, 2001). In this study (see 
Table 1), the nucleotide Tamura-Nei distances between each of four haplotypes (I to IV) of N. 
caudatus in Korea and the haplotype (V) of N. caudatus in China varied from 0.0650 to 0.0803. 
It indicates that the level of cytochrome b sequence diversity in N. caudatus is high in Bovidae. 

The habitat fragmentation of wild animals is caused by human industrialization, and the reduction 
of genetic diversity or extinction occurred rapidly within isolated areas (Birky, 1991; Lande, 1993). 
The last decade has seen the beginning of a field that applies the principles and methods of 
population genetics to species conservation, and DNA sequence analysis is now routine (O'Brien, 
1994; Cassens, 2000). Barratt et al. (1999) found one or two mitochondrial DNA control region 
alleles within 11 of 15 populations of red squirrel (Sciurus vulgaris), and they supposed that many 


18 Korean J. Syst. Zool. 18(1), April 2002 


0.01 
57 [EN 
96 

Iv 

N 

S 

76 y 

83 [ 

vi 


Fig. 3. Phylogenetic trees with six haplotypes (276 bp) in two species of Naemorhedus. Trees were 
constructed by neighbor-joining (A) and maximum parsimony (B) with 500 bootstrapped replications. Tamura- 
Nei distances presented in table 1 were used. Haplotypes are labeled as in Fig. 1, and the sequence (S) of Saiga 


tatarica was used as an outgroup. 


population experienced demographic contraction in the last hundred years. Sokal and Oden (1978) 
noted that a continuous increase of genetic distance with geographic distance can be expected 
under isolation by distance. 

Grey goral (N. caudatus) distributes in extreme southeastern Siberia to Burma and western 
Thailand, but International Union for Conservation Nature (IUCN) classifies it as vulnerable because 
of poaching and destruction of forest habitat (Nowak, 1999). In this study with N. caudatus in 
Korea and China (see Table 1 and Fig. 2), the Tamura-Nei distances among three haplotypes at 
Yanggu were 0.0151 to 0.0185, and it suggested that the genetic diversity of Yanggu population 
was decreased in low level. And the distances between each of three haplotypes at Yanggu and 
one haplotype at Samcheok were 0.0343 to 0.0489, and the distances between N. caudatus 
populations of Korea and China were 0.0650 to 0.0803, indicating that genetic distances 
appeared to be increased with geographic distances by habitat isolation. Moreover, it was also 
revealed that the genetic diversity of N. caudatus in Korea was decreased, and various 
conservation plans to mitigate the reduction of the diversity have to be in immediate action. 

Mitochondrial DNA is a valuable genetic marker for studies of evolutionary relationships of 
species (Sunnucks, 2000). N. caudatus was a subspecies of N. goral (Tate, 1947), but Groves and 
Grubb (1985) considered it to be a distinct species. In this study with mtDNA cytochrome b gene of 
two species in Naemorhedus (Table 1 and Fig. 3), it is concluded that N. goral was not distinct 
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from N. caudatus. It suggests that they are conspecific, as noted by Tate (1947), but further 
analyses with additional speciemens of two species are necessary. Cronin (1991) noted that 
relationships based on mtDNA might not be concordant with true species phylogenies in Cervidae. 
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